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ABSTRACT 

In a recent analysis of the radio emission from the planetary nebula A30, Dgani, Evans & 
White (1998) claim that the emission, located in the inner region, is probably dominated by 
nonthermal emission. We propose a model to explain this. We assume that the fast wind, blown 
by the central star of A30 carries a very weak magnetic field. The interaction of this wind with a 
cluster of dense condensations traps the magnetic field lines for a long time and stretches them, 
leading to a strong magnetic field. If relativistic particles are formed as the fast wind is shocked, 
then the enhanced magnetic field will result in nonthermal radio emission. The typical nonthermal 
radio flux at 1 GHz can be up to several milli-Jansky. In order to detect the nonthermal emission, 
the emitting region should be spatially resolved from the main optical nebula. We list other 
planetary nebulae which may possess nonthermal radio emission. 

Subject heading: ISM: magnetic fields — ISM: kinematics and dynamics — Radio continuum: 
ISM — Planetary nebulae: general — Planetary nebulae: individual (A30) 
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1. INTRODUCTION 

The possibility of magnetic fields carried by the fast winds in planetary nebulae (PNs; 
Chevalier and Luo 1994; hereafter CL; Chevalier 1995; Garcia-Segura 1997) raises the question 
of whether the interaction between this wind and the previously ejected slow wind gives rise to 
particle acceleration, and hence nonthermal emission from relativistic electrons (Chevalier 1995). 
Chevalier notes that the major problem with detecting the nonthermal emission is the strong 
thermal free-free emission emanating from the dense shell. In many PNs, however, knots and tails 
are present relatively close to the central star, where both the ram pressure of the fast wind and 
the magnetic field are higher, increasing the nonthermal emission. In the hydrogen deficient PNs 
A30 and A78, for example, knots with tails appear in a disk structure close to the central star 
(Borkowski, Harrington, Sz Tsvetanov 1995). These nebulae are believed to have had a helium 
flash late in their evolution, and are known to possess fast winds, v w ~ 4000 km s _1 (Leuenhagen, 
Koesterke & Hamann 1993 ). The knots are believed to be photo evaporating globules (Dyson 
1968; Kahn 1969; Borkowski et al. 1995). The HST images of the knots in the central region of 
A30 (Borkowski et al. 1995) show that they are being disintegrated and accelerated by the dilute 
shocked fast wind from the central star, giving rise to Rayleigh- Taylor instability. In addition, 
hydrogen deficient nebulae have lower thermal emission (see discussion in Dgani, Evans & White 
1998), allowing the nonthermal component to be detected. 

Jones, Kang & Tregillis (1994) have simulated the evolution of dense gas clouds moving 
supersonically through a magnetized low density medium. They took a passive magnetic field and 
found that the synchrotron emission from relativistic electrons increase rapidly when the cloud 
begins to fragment under Rayleigh- Taylor instability. The synchrotron emission is dominated by 
field enhancement through stretching in the shear turbulent layers formed between and behind the 
fragments. Jones, Ryu & Tregillis (1996) performed full two dimensional MHD calculations and 
found that when the ambient field is transverse to the flow fast amplification of the field occurs 
also at the nose of the cloud, and not only in the turbulent regions. 

Dgani et al. (1998) present a radio map of the central region of A30. The map shows an 
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extended source of about 0.15 mJy stretching from the central star to a point 3" to the west. 
Based on a comparison to the optical line emission from that region, they claim that the source is 
probably nonthermal. Motivated by this recent finding, we apply the theoretical results of Jones et 
al. (1994; 1996) to the interaction of a magnetic fast wind with knots in PNs. The amplification of 
the magnetic field is calculated in §2, while in §3 we estimate the nonthermal radio emission from 
the interaction region. In §4 we summarize, and list PNs, other than A30 (Dgani et al. 1998), 
which we expect to have nonthermal radio emission. 



2. INTERACTION OF DENSE KNOTS WITH A MAGNETIZED WIND 

Like CL, we examine a slowly rotating central star, for which the tangential component of 
the magnetic field dominates far from the star. For a field line that starts radially from the stellar 
equator, the downstream tangential component of the field at r > r s , where r is the distance from 
the central star and r s its radius, is given by 

B = B S V -^^, (1) 
v w r 

where t> ro t is the equatorial velocity of the rotating star, v w is the fast wind terminal velocity, and 
B s is the (radial) magnetic field on the stellar surface. Since both the magnetic pressure B 2 /Sir 
and the ram pressure pv^ drop as r~ 2 before hitting the shock wave, it is useful to define the 
parameter (CL): 

B 2 B]r[ fvruV 



(2) 

4irpvl M Vw V v w J 

where M is the mass loss rate to the fast wind and p = M(^itr 2 v w )~ l is the density in the wind, a 
is the ratio of field energy density to kinetic energy density in the wind. CL find that the magnetic 
field shapes the nebula for a ;> 10~ 4 . From equation (2) we can write the magnetic field in the 
equatorial plane in the free wind as 

B = 3 5 x 10- 7 1 ( ) ^ ( Vw \ 1 ( T \ 1 C (S\ 

UO-V I 5 x 1O~ 8 M yr- 1 ) UoOOfem s~V V i0 16 cmy ^ W 
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The mass loss rate is scaled according to the value Borkowski et al. (1995) used for A30. 

Consider first a single globule with radius R g . A supersonic flow incident on it will form a 
bow shock at its nose, a turbulent shear layer around its sides and a turbulent wake behind it. 
The magnetic field is expected to be enhanced in these regions because of line stretching (e.g., 
Jones et al. 1994, 1996). As shown in the extensive numerical simulations of Jones et al. (1996), 
when the magnetic field lines are transverse to the flow the magnetic pressure on the nose of the 
clump reaches equipartition with the ram pressure, even if the initial magnetic field is quite weak. 
Using their analysis for the numerical results, we can estimate the growth rate of the magnetic 
field by the following consideration (their equations 8-10) 

Jones et al. (1996, eq. 8) start with the ideal MHD magnetic induction equation 



u, p and B are the velocity, the density and the magnetic field in the flow respectively. Taking 
(B ■ V)u ~ Bu/Rg, equation (4) becomes 



The last equation gives exponential growth to the magnetic field as long as the magnetic field 
lines stay near the stagnation point. Jones et al. (1996) estimate that at the stagnation point the 
magnetic pressure will rise exponentially until the bullet is crushed or until the magnetic pressure 
equals the ram pressure. However, since they performed 2D calculations, they neglected the sliding 
of the magnetic field line around the bullet. Although they claim irregularities in the 3D spherical 
bullet will capture the field lines and prevent sliding, we believe that for a flow with very large 
(3 = -Pth/-PB> the field lines will be carried around the globule in a time <^ R g /u, much shorter than 
the bullet's crushing time. Here P t h and Pg are thermal and magnetic pressures, respectively. 

Other sites of magnetic field enhancement are the turbulent mixing layer at the sides of the 
globule, and its turbulent wake. These are more promising sites, mainly because the clumps 
(globules) appear in groups, and turbulent flows dominate the region. 



dln(B/p) _ B-[(B- V)u] 
It ~ B~ 2 



(4) 



dln(B I p) u 
It ~ ~R g ' 



(5) 
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In here, we consider the flow around and inside an ensemble of clumps, occupying a region of 
size i?ens- The typical distance between the clumps is of the order of their size R g . This case has 
not been studied numerically, but we can apply an analytical approach to estimate the magnetic 
field enhancement. A turbulent flow will form between the globules, where the magnetic field can 
be amplified. There, the change of the wind velocity u along a field line depends on the size of the 
mixing layer I ~ R g . That is (B ■ V)u ~ Bu/R g . The time a field line spends in the turbulent 
region is the time it takes the flow to pass the group: R ens /u. Substituting these values in the 
induction equation (4), and integrating over the time the field lines spend inside the group of 
globules, gives the relative amplification 

B/p ~ exp^p?-. (6) 

Rg 

For R g ~ 0.1i? ens , for example, there is a large amplification factor ~ 10 4 . Therefore, the magnetic 
field will reach equipartition even if its initial value is ~ 10 -4 that of equipartition, i.e., a ~ 10~ 8 
(eq. 2). In the model of CL and Garcia-Segura (1997) the post shock magnetic field grows as 
B ~ r, which gives a factor of ~ 10 — 100 for the magnetic field and ~ 100 — 10 4 for the magnetic 
pressure. Since CL assume that the magnetic field shapes the PN globally, they have to assume 
that the initial magnetic pressure is ;> 10~ 4 /w 2 , where pv 2 is the fast wind's ram pressure. This 
means a > 10 -4 . In the flow pattern assumed here, the localized turbulence results in a much 
larger increase in the magnetic field. Therefore, to achieve equipartition between the magnetic 
pressure and the post shock thermal pressure, a value of a ;> 10~ 8 is sufficient. 



3. NONTHERMAL RADIO EMISSION 

The upper limit for the enhanced magnetic field is the equipartition value, found from the 
relation Pg = pv^, 

B max ~ 5 x 10~ 3 ^- - w _ 8Mq ( 4000 km 8 -i) (iq16^) G ' (7) 

where d is the distance between the central star and the group of globules. If the relativistic 
electrons are produced in the bow shock and then carried downstream to the magnetized region, 
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they can emit synchrotron radiation. We assume equipartition between the magnetic field and 
relativistic particles (electrons and nuclei). The radio luminosity between v\ and v 2 can then 
be estimated if the spectral index a, (i.e., the flux F u oc v~ a ), of the radiation is known (cf. 
Pacholczyk 1970 Page 171, Eq. 7.14). The radio luminosity is given by 

B 7 / 2 $ R 3 

''ens (q\ 

~ 4.5 (1 + fe) c 12 (a,v 1 ,v 2 y {) 
where $ is the fraction of the ensemble volume occupied by magnetic field and relativistic particles, 
k is the ratio of the heavy particle energy to the electron energy, and Ci2(a, V\, v 2 ) is tabulated in 
appendix 2 of Pacholczyk. 

Taking v x = 10 7 Hz and v 2 = 10 11 Hz, a = 0.75, for which c 12 = 5 x 10 7 and k = 100, 
(Chevalier 1977; 1995) we obtain 

L ~ 4.9 x 10 28 3> ( Ren ° ) 3 ( ^-_y /2 er g S -i. (9) 

V5xl0 15 cmy V5xlO- 3 G/ & v; 

The radio flux at 1GHz can now be estimated: 

Fv9 ~ 3.5$ ( Ren ° ) 3 ( V2 ( —) ~ 2 mJy (10) 

V5xl0 15 cmy \5xW- 3 Gj \kpcj y v ; 

where D is the distance to the source, and lmJy = 10~ 29 W vaT 2 Hz -1 . A few mJy is much less 
than the observed fluxes from the shells of most catalogued PNs (Zijlstra, Pottasch, & Bignell 
1989). However, the emission we predict is far from the PN shell, near the central star, and will 
not be overwhelmed by thermal emission. Chevalier (1995) discussed the case in which nonthermal 
electrons are produced in the fast wind termination shock, and then carried downstream to the 
magnetized region near the shell. In this case the nonthermal emission region neighbors the 
thermal emission region and cannot be separated observationally. Therefore, Chevalier concluded 
that a search for radio synchrotron emission from planetary nebulae would involve a frequency 
lower than 1GHz. 

Dgani et al. (1998) show an 8 GHz radio map of the hydrogen deficient PN A30. They 
present a source of ~ 0.1 mJy which extends to about 3" west of the central star. A30 possesses 
a large shell (120") visible in both H Q and OIII, and a hydrogen deficient central region visible 
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in OIII but not in H a (Jacoby 1979). The radio emission emanates from the central region. The 
source extends to the vicinity of the main blobs of the HST image (Borkowski et al. 1995), but it 
cannot be identified with any of them. This is compatible with our model where the radio emission 
results from magnetic amplification in the turbulent region formed by decelerating clumps. 

A somewhat parallel situation occurs in young SNRs. In the SNR Cassiopea A (Cas A), the 
radio knots as well as the optical knots appear to be connected to the clumpy component of the 
ejecta, however, they are not identical. Anderson et al. (1994) analyze the observations of optical 
and radio emission knots in the SNR Cas A. They claim that while the clumps are almost certainly 
sites of particle acceleration, when there is substantial magnetic field growth, it will enhance the 
synchrotron emission. Complex flow patterns bring in particles that were accelerated in different 
regions. In young SNRs particles are produced in many acceleration sites: the blast wave shock, 
the reverse shock and many internal interactions shocks. Thus the enhanced emission in the radio 
knots is probably a result of field amplification. The numerical simulations of Jun & Jones (1998) 
and of Jun, Jones & Norman (1996) support this conclusion. 



4. DISCUSSION AND SUMMARY 

We studied the interaction of a weakly magnetized fast wind, of velocity v w > 1,000 kms -1 , 
blown by the central stars of PNs with high density condensations, and we propose that it leads 
to the amplification of the magnetic field. When this interaction is with a cluster of clumps, the 
magnetic field lines are trapped for a long time in the interaction region. This process is equivalent 
to the numerical simulations of Jones et al. (1996), in which the 2D nature of the flow did not 
allow the magnetic field lines to slide around the condensation. Jones et al. (1996) found that the 
magnetic field can reach equipartition with the wind's ram pressure. 

The enhanced magnetic field will have several effects: (1) It will reduce the mixing 
(entrainment) of the clump material with the fast wind. (2) It will reduce the heat conduction 
between the cool clump and the hot shocked fast wind material. (3) The magnetic field lines will 
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reconnect among themselves. This process may lead to heating and further X-ray emission. (4) It 
may result in nonthermal radio emission. The first three processes will affect the X-ray emission 
from the hot material. In the current study we did not examine these processes. We explored 
here, and concluded that the expected nonthermal radio emission can be detected in PNs with fast 
winds and clumps close to the central star. Such a PN is A30 (Dgani et al. 1998), claimed also 
to have extended X-ray emission partially overlap with the location of the optical condensation 
(Chu, Chang & Conway 1997). 

The required initial magnetic field for detectable nonthermal radio emission is ~ 2 orders 
of magnitude lower than the field intensity required to influence the shape of PNs, according to 
the model of CL and Garcia-Segura (1997). The stronger magnetic field required by CL and 
Garcia-Segura (1997) requires much faster rotation, as shown by equation (2) with a > 1CP 4 . This 
is very problematic for their model, since the central stars of PNs are expected to rotate slowly 
because of mass loss. This is true even if their AGB progenitors are spun-up by companions on 
the AGB (Soker 1998). The field in the model of CL and Garcia-Segura (1997) is globally large. 
The amplification of the magnetic field through the interaction of the magnetized wind with a 
cluster of condensations, on the other hand, enhances the field locally, but not globally. Therefore, 
the detection of strong magnetic fields in specific locations in PNs does not necessarily mean that 
the magnetic field is globally strong. 

The process proposed in this study predicts detectable radio synchrotron emission from PNs 
having fast winds and high mass loss rates, both of which determine the magnetic field intensity 
(eq. 7), and optical condensations. In order to detect this synchrotron emission we need to resolve 
the central region of the nebula from the shell. This requires high spatial resolution observations. 

To list several PNs, in addition to A30 and the other hydrogen deficient PNs (see discussion 
in Dgani et al. 1998) , for which we expect nonthermal radio emission we select PNs having fast 
winds from Patriarchi & Perinotto (1991) and check their optical images for condensations in the 
inner region. We find that from the list of Patriarchi & Perinotto (1991) only NGC 2371 and NGC 
7094 show no structure in the very inner regions. Therefore, in general all PNs with fast winds 
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blown by their central stars should be searched for nonthermal emission, in particular, those with 
strong condensations in their central regions: A78, NGC 5189 and NGC 6210. 
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